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Abstract 

In this paper we discuss analogue computers based on quantum optical systems ac- 
celerating dynamic programming for some computational problems. These computers, at 
least in principle, can be realized by actually existing devices. We estimate an accelera- 
tion in resolving of some NP-hard problems that can be obtained in such a way versus 
deterministic computers. 
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1 Introduction 



Quantum computing is intensive developing intersection of physics and informatics since 
(Feynman 1982; Shor 1994; Grover 1997; Preskill 1998). But schemes of universal quantum 
computers working by qubits meet formidable difficulties in their realization as real physical 
devices (Dyakonov 2003). In this note we develop another approach based on analogue re- 
alization of quantum computing. We concentrate our attention on some particular NP-hard 
problems. For a given problem we can try to construct a quantum analogue machine that 
resolves this problem by means of time evolution. 

Recall (Garey and Johnson 1979) that, roughly speaking, a problem with instances I 
lies in the class NP if there is a polynomial time _F*( |_T| ) algorithm checking a solution (if 
this solution is provided) where |7| is the size of the input I. The problem is NP-hard if 
P = NP whenever this problem belongs to P. The list of NP-complete problems contains 
many important problems of number theory, graph theory, logic etc. Some of the most known 
NP-complete problems are travelling salesman problem, satisfiability, knapsack, matrix cover 
etc. (see a list in (Garey and Johnson 1979)). Let us observe that all the NP-complete 
problems, are, in a sense, polynomially equivalent. This means that there are polynomial 
time reductions between different problems. 

An analogue quantum machine or resolving of combinatorial search NP-hard problems 
has been proposed by (Farhi et al. 2000). Let us consider 3-SAT problem in which we deal 
with a family of N clauses of the form (B^ \J B; l2 \j -< £>j 3 ) or (B^ \J -> B{ 2 \J -> J5j 3 ) or 
(Bi x V Bi 2 V Bi 3 ) etc., where Bj are boolean variables taking values True or False. The 
problem is whether one can assign such values to all the variables that all the clauses will 
become true. 

The approach of (Farhi et al. 2000) uses quantum adiabatic theorem. The authors 
describe a formal Hamiltonian such that its ground state gives us a solution of 3-SAT problem. 
The bound of the actual speed-up yielded by this algorithm depends on the spectral gap 
between the ground state and the first excited state. If this gap is exponentially small the 
time of solving is exponentially large. It is a very difficult problem to estimate the value of 
this gap and thus it is not easy to show, in general, that such an algorithm really gives a 
quantum acceleration. (Another NP-complete problem that can be naturally associated with 
a Hamiltonian is the problem of minimization of energy of spin glass with a large number of 
spins (see (Garey and Johnson 1979, p. 282.)) 

The second difficulty with such approach is that it is not obvious to realize physically 
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qubits and Hamiltonian with very nonlocal nontrivial interaction between qubits. It seems 
that difficulties in physical realization of Feynman quantum computers or analogue computers 
from (Farhi et al. 2000) are significant and some authors even believe (see, for example, 
(Dyakonov 2003)) that such computers are physically non-realizable. 

The approach of (Farhi et al. 2000) uses a specific structure of 3-SAT problem. If we 
use some universal approach that does not take into account a specific form of the problem, 
we can expect a quantum acceleration in time iV 1 / 2 (Grover 1997), i.e., if a deterministic 
computer makes a search in N steps then quantum universal computers will make the same 
search in N 1 / 2 steps. Further, it was shown that acceleration, with a universal approach can 
be achieved by at most iV 1 / 6 Beals et al. 1998; Preskill 1998. Notice that this quantum 
N 1 / 2 acceleration is less than an acceleration that can be obtained by special deterministic 
and probabilistic algorithms with respect to trivial exhaustive search (Dantsin et al. 2003). 
Once more quantum machine based on quantum optics phenomena was discussed in (Kazakov 
2003). The earlier analogue computational machines were discussed in (Matiyasevich 1987; 
Blass 1989). 

In this paper we describe a quantum machine which, as we believe, can be realized phys- 
ically and which may accelerate solving of some NP-complete problems. Our machine uses 
different linear and nonlinear quantum optical devices that can be constructed actually and, 
thus, our computer, at least in principle, can be practically realized. An acceleration that 
can be obtained heavily depends on different characteristics of our devices. We estimate this 
acceleration vs. deterministic computers using the characteristics of actually existing devices. 

2 Statement of problems. Known results. 

We consider the following problems which are NP-hard (Papadimotriou, Steiglitz 1982; Garey, 
Johnson 1979). 

1 Boolean knapsack, variant 1 

Given positive integers Cj, j = 1, ...n and K, is there a subset S of {0, 1, ...,n} such that 
= JT7 

In other words, whether there exist n boolean values si € {0, 1} such that Ya=i CiSi = ^ 
Here the size |/| of the instance {n, c±, C2, Cn, K} is the number of bits needed for 

binary representations of all the integers ci,K. We can suppose, without loss of generality 

that q < K. Thus, size |/| can be estimated roughly as 0(n log K). 

2 Boolean knapsack, variant 2 
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Given integers Cj, j = l,...n and B_,B + , whether there exist n boolean values Sj G 
{0,1} such that Ya=1 °i s i nes i n the interval (B-,B + )? Here the instance size is roughly 
0(n log B+). 

3 Optimization boolean knapsack 

Given integers Cj and iuj , j = 1, ...n and the number B+, maximize the cost 

n 

(1) 

1=1 

defined by n boolean variables s\ £ {0, 1} under condition that 

n 

Y^c iSi <B + . (2) 

i=\ 

There is an important difference between the problems 1, 2 and 3. The output in 1,2 is 
"YES" or "NO". The output of 3 is a number, and one could try to approximate it. 

Let us remind now some known results about 1, 2 and 3. 

a All the problems 1-3 can be resolved by exhaustive search in 2™ steps. 

b If 2 n > K then the problem 1 can be resolved by a more effective method, namely, 
by dynamic programming, in 0{nK) steps. This method can be described briefly as follows 
(for details, see Papadimotriou and Steiglitz 1982). The algorithm produces consecutively 
Eo, Si, E n such that E,- is the set of all possible subsums of c±, Cj At the first step we 
set Eo = {0}. At j + 1-th step, we set Ej+i = (Ej U (Ej + Cj + i)) n {0, K}. The problem 
has a solution if at the last n-th step K £ E n . 

In a similar way, we can resolve problem 2, it takes 0{nB + ) steps, and problem 3, it takes 
0(n min(i?_|_, R op t)) steps, where R op t denotes the maximum cost (0). 

c Suppose we solve an approximative problem 3, namely, we seek for a number close 
to a maximal cost. This means that we give up accuracy in exchange for efficiency of our 
algorithm. To this end we can apply a truncation procedure removing the last t digits from 
the binary representations of wi and Cj . Let w m be the largest W{ . Such a procedure leads to 
a cost Rappr satifying the estimate (Papadimotriou and Steiglitz 1982; Ibarra and Kim 1975) 

Ropt Rappr ^2 

p ' e = 7T~ ( 3 ) 

I^opt uJ m 

The approximating algorithm runs in time 0(n 4 e _1 ) (Papadimotriou and Steiglitz 1982; 
Ibarra and Kim 1975). Thus this problem has an approximative solution that can be found in 
a polynomial number of steps. Notice however, that there are many NP-complete problems 
that cannot be approximated in such a way, for example, travelling salesman problem (Garey 
and Johnson 1979; Papadimotriou and Steiglitz 1982). 
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3 Description of the quantum machines 



Consider n + 1 points xq < x\ < X2, ■■■ < x n = xt located along x— axis in (x,y) -plane. At 
the first point we set a laser, which generates a narrow beam. The diameter of this beam will 
be denoted by d&, the wave length of the laser radiation will be denoted by A. Typical values 
of these parameters are df, ~ 2 • 10 -3 m, A ~ 5 • 10~ 7 m. 

Here we describe an analogue quantum optical device (QOD) for the knapsack problems 
1,2. Its possible scheme is presented on fig 1. An input laser beam is splitted by 50% mirror 
Al to two separate beams running two different trajectories. The beam 1 then is shifted by 
a plane optical plate on the value c%k in (vertical) z-direction, which is perpendicular to the 
plane of our picture (here k means the minimal shift). Then the beams are on the second 
mirror Bl and united result (which is a combination of shifted and non-shifted beams) goes 
through amplifier CI. We suppose that this amplifier has the characteristics presented on fig. 
2. After passage of m mirrors the propagating light will contain beams shifted in z-direction at 
all possible distances YliLi c i s i K i where the values G {0, 1} depend on the trajectory of the 
beam. At the final point we measure intensity of outcoming light by a charge-coupled device 
(CCD). A CCD camera uses a small rectangular piece of silicon to receive incoming light. 
This silicon wafer is a solid state electronic component which has been micro-manufactured 
and segmented into an array of individual light-sensitive cells called "pixels". The pixel of 
the most common CCD has size only about 6 P ~ 10~ 7 m ~ 0.2A, it measures the intensity of 
light independently from other pixels. Note, that namely 50% mirrors are genuine quantum 
components in our device (Feynman 1985). 

We denote by Rm the z-size of the separating mirrors (and amplifiers). The last parameter 
important for estimating of quality of our machine is an angular divergence a of the laser 
beam. A usual laser produces a beam having the angular divergence of order a ~ d&/A ~ 
4- 10~ 5 . 

Our machine is completely defined thus by the following parameters: 

^ j Rm > 3) ^min 5 ^max j A, df, , Oi , fC, Cj . 

For the problem 3, we use the following ideas. In order to realize sums (0 |2J) we have to 
modify scheme presented on fig 1. Namely, together with z-shifts we use the horizontal beam 
shift in an orthogonal direction y. Notice that a beam shifted in y-direction contains this shift 
propagating on different trajectories L\ and L2 (see fig. 3). It means that after successive z- 
and y-shifts, which the initial beam gets propagating from xq to x n , one obtains the set of 
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beams whose z- and y-shifts are different sums 

n n 
i=l i=l 

in accordance with trajectories. 

At the end of the device we set CCD which measure the plane distribution of the outcoming 
light with z-coordinates z < B + — 5 P , where 8 P is as it was mentioned above the typical pixel 
size. 

4 Solving the dynamic programming by the quantum optical 
device 

Consider possible laser beam trajectories as a result of n reflections on our mirrors. After 
the first series of reflections we obtain two beams running along lines y = h\ and y = since 
each photon randomly chooses a way either along the first trajectories or the second one. So, 
the mirror system serves as a gate choosing the photon way. Remark that the mirror can be 
considered simultaneously as a quantum and a classical device. In a certain sense, our gates 
can be thus named as semiclassical quantum gates. 

If we set a CCD that registrates photons at a point located after xi, we could registrate 
two localized z-separated gaussian beams at y = and y = H\. 

So, in an ideal situation, QOD works as follows. Consider the problem 2. We set Hi = nci, 
i = 1,2, where Cj,n are given. We take now Y m i n = kB-, Y max = kB + . Now turning 
on an input laser beam, we check, does the interval (Y m i n ,Y max ) contain the center of the 
any laser beam (respectively, answer either "YES" or " NO"). We suppose that CCD solves 
the problem with an error less than 1/3. 

It is not difficult to see that this method of solution can be considered as a physical 
realization of dynamic programming method from Section 2. 

The problem 3 can be resolved approximatively, i.e., the global maximum can be estimated 
within some precision. However, to this end, we must modify slightly our consideration (see 
below) . 

Consider now diffraction effects. They can create an obstacle since one could registrate a 
photon not really passing correct gates as a photon that resolves our problem. In accordance 
with theory of gaussian beams (Svelto 1982) each gaussian beam has a finite transversal size 
and an angle of divergence a. These values are connected by 

a ~ ——. (4) 

d b 
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After n reflections, the size of an output beam becomes 

d final = nL sin a + d b , (5) 

where L means distance between mirrors. Therefore, to get minimal value of dfi na i we set 

db = VnLX. (6) 

Further, in order to separate on registering CCD the broadening gaussian beams (whose 
amplitude can differ up to factor 2 in accordance with fig.2), we have to restrict the minimal 
distance k between beams by 

k > df ina i = 2(4- (7) 
In addition to (4.3) we have a geometric restriction 

Rm > ku + d . (8) 

This inequality means that the shifts of no beams jump out any mirrors. 

One more problem is a possibility that, at some step i = Tq many different beams arrive 
at the same point. It is not a mathematical difficulty (since this means the existence of many 
solutions) but it is a physical difficulty since it can lead to a high energy concentration and 
destruction of the mirrors. To avoid this effect we use active media that saturates the photon 
density (see fig. 2). This allows us to restrict the photon density at each point of the mirror 
by some constant. 

Using of the active media leads however to new difficulty which can give us a lost of 
solutions. Namely, the phases of the photons of each beam become slightly different after 
passing the active media. Thus if the problem 2 has many solutions, CCD could not registrate 
any photons as a result of interference. This effect is possible if a phase shifts are significant. 
To overcome this difficulty we use special auto phase control devices known in optics. At the 
points z\,Z2, ■■■Zn (see fig. 1) we places the phase-adjusting devices which measure phase of 
all beams and compare it with phase of a reference laser beam of the same frequency. Such 
adjustment gives the possibility to align phases of mixing beams. 

5 Estimation of QOD performance for boolean knapsack 

Let us estimate now what we can do using such a QOD and compare this machine performance 
with deterministic computer performance. 
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5.1 General estimates 



There are three important parameters: a) the implementation cost; b) the energy cost of 
solving; c) the running time when we solve M times the same problem with different inputs. 

We compare here the dynamic programming for knapsack problems 1 , 2 by a deterministic 
computer and by QOD. For the sake of unifying the notations we assume, that B + < K. 

For the QOD the implementation cost CI quan t can be estimated roughly as 

CI quant = 0(Kn) 

since the mirror and amplifier sizes are of order K , the number of the mirrors and the cuvettes 
is n. The amplifier size is of order Kn, one can expect that energetical cost CE q i per one 
instance of the problem also is 

CE qi = 0(Kn(n + K)), 

where factor n + K arises due to greatest length of the photon trajectory. The complete 
energy involved in M calculations is then 

CE quant = 0(Kn{n + K)M). 

The implementation time has the order Kn and the resolving time 

Time quant leC\(n + K)M, 

where C\ is inverse proportional to the light speed and thus this coefficient is rather small. 
For the deterministic computer we have 

CE det = O(KnM), Time det = O(KnM) 

Thus, QOD wins in time (taking into account both implementation and running times) vs. 
deterministic machine, while the order of consumed energy is higher than for a deterministic 
one. 

5.2 Energetic and temporal costs 

Our results include energetic and temporal costs of calculations. Let us consider a model, 
where this correspondence exhibits in a more explicit form. Let several polyhedra be placed 
in R 3 and we have to calculate the length of the shortest path from the point A to the 
point B avoiding these polyhedra (treated as obstacles) (see fig.4). It is well known that this 
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problem is NP-hard (Canny and Reif 1987). Moreover, even determining first 0(N 1 / 2 ) bits 
of the length of the shortest path is NP-hard, where N denotes the bits size of description of 
the instance of the problem (Canny and Reif 1987). We can realize this calculation by the 
following analogue physical device. Let us place in the point A a source of radiation, in the 
point B - detector and polyhedral obstacles have mirror boundaries. If we turn the source 
A at initial moment t = and measure the time t = t\ when our detector catches a first 
photon. Let the power of source be J and frequency of radiation be v, then the number of 
created photons is J /his, where h is Planck's constant (Karlov 1992). Each photon moves 
in its path and one can describe our analogue calculation as a "parallel computation" of the 
shortest path by different photons in which each photon plays the role of a processor. 

In general, this first photon caught by the detector could not belong to the front wave 
(which corresponds to the shortest paths) due to a weakening of the wave and possible presence 
of the more strong waves. Therefore, increasing the intensity of the source allows to augment 
the probability to detect a photon of the wave front. So, we have to either repeat this 
experiment, with corresponding increasing of temporal cost, or increase the intensity of the 
source with increasing of energetic cost of calculation. 

On the other hand, considered analogue devices have brought us to a conjecture that a 
tradeoff of the following kind 

E • T > Z, 

holds, where E and T are energetical and temporal costs of this calculation respectively, and 
Z depends on the problem. This means (at some extent) that there could be an "exchange" 
between the energy cost and the temporal cost needed to solve a computational task. 

5.3 Estimates for real parameters 

For the rate of a deterministic computer, for a single processor (we do not use a parallelism) we 
take the value V = 10 10 operations / sec. Suppose 2™ > B + , then the dynamic programming is 
more effective than the exhaustive search. Then a deterministic processor solves the problem 
2 within time T c 

T c « nB + V~ l (9) 
Consider our optical processor. We see that B + must be subject to 

B + < R m /k. (10) 
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Thus our possibilities are restricted by the mirror size. Then the problem can be resolved 
within time 

T q pa (nh + R M )/c+ nT atom , (11) 

where the second term describes the time when photons spend in the amplifier (which is 
at most n times the relaxation time T a t om for active atoms, T a t om is about 10~ 8 sec). The 
preprocessing (implementation) time for the optical machine is 

Tq^pre ~ C2K7I, 

where C2 is a constant. 

It is interesting now to see what we could obtain by using really existing devices. As an 
example, we take the following typical parameter values: 

R M = n L = 10m, d b = 2 • 10~ 3 m, A = 5 • lCT 7 m. 

Then k = 5 ■ 10~ 3 m and B + < 2 • 10 2 . For n > 30 the next inequality holds, 2" >> nB + , 
and dynamic programming is more effective than the exhaustive search. 

Minimal admissible difference Y max — Y m in can be estimated 5 • 10~ 3 m. The 

deterministic time (5.1) is then T^et ~ 30 • 2 • 10 2 • 10~ 10 ~ 10 _6 sec. In accordance with (|llj) 
T q pa 30/(3 • 10 8 ) + 30 • lO" 8 pa 10~ 6 sec. 

So, our conclusions for problem 1, 2 are the following. The performance of QOD is 
restricted by mirror sizes and beam diameter. The problem 2 can be solved by QOD and 
this effectiveness is essentially more vs. a deterministic processor if we repeat the same 
computation with different inputs many times. For realistic values of parameters, considered 
above, we can handle the case n < 60, and then the speeds of the QOD and deterministic 
processor have the same order. 

Let us discuss now situation for problem 3. In the next section we describe first modifi- 
cations needed for resolving this problem. We will see that in this case the speed-up is much 
better. 

6 Estimate of machine performance for optimization quantum 
knapsack 

6.1 General estimates 

We compare here the dynamic programming for knapsack problem 3 by a deterministic com- 
puter and by QOD using the parameters CI, CE, Time described in Section 5.1. 
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For the quantum optical device the implementation (preprocessing) cost CI quan t can be 
estimated roughly as 

CI quant = 0(K 2 n), 

since the the mirror size is now of order K 2 , the number of the mirrors is n. The amplifier 
takes volume of order K 2 n, one can expect that the energetical cost CE quan t is 

CE quant = 0{K 2 n(n + K)M), 

where now K majorates sums and (jSJ). However, if only numbers Wi are large and c% ~ 1, 
then the mirror area becomes again O(K) as in Section 5, and we have the same estimates 
as in subsection 5.1. 

One can estimate the resolving time 

Time quant = 0({n + K)M). 

For the deterministic computer we can (if we solve M times the problem with different 
inputs) estimate energy similarly to 5.1 (problems 1,2) 

CE det = 0(KnM), Time det = 0(KnM). 

Thus, if the number of inputs M » K then QOD has a chance to win, in time, with respect 
to a deterministic processor. The consumed energy of the QOD can be estimated by 

CE quant = 0{K 2 n(n + K)M) 

(cf. subsection 5.1 above) whereas the energy consumed by a deterministic machine can 
be bounded again by O(KnM) (cf. the discussion on the tradeoff of energy and time in 
subsection 5.1). 

Let us consider now approximating solutions (see Section 2). Recall that the processing 
time of the deterministic processor in order to obtain an approximative solution, within 
precision e, is 

Time detjappr = ©(MtA" 1 ). 

To compare this performance with effectiveness of our QOD, let us note that the pixel size 
is 5 p and thus the size of the mirrors in QOD solving the problem within relative precision e 
should be 5 p /e = Kk. Then the implementation cost will be 

CIquant = 0((5 p / eK) 2 n) , 
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resolving time can be estimated as Time quan t = 0((n + 5 p /en)M) The choice S p is restricted 
by a small number of possible values, thus, Time quan t + Cl quan t grows in n more slowly than 
Time de t,a PP r for e _1 = o(n 3 M). 

The consumed energy of QOD for approximating problem can be estimated by 

CE quant:appr = 0((5 p /eK) 2 n(n + S p /eK)M) 

which is less than the consumed energy for a deterministic machine (which is of the order 
n 4 e _1 M) when e" 1 = o{(nn/ 5 p fl 2 ). 

6.2 Estimates for real device parameters 

Recall that, to resolve the problem 3, we have modified QOD introducing additional parallel 
plates performing beam shifts in the direction y. So, we can suppose that i-th gate makes 
the shift to Hi in the direction z and to H[ to the y axis. As above, we set Hi = KCi and 
H[ = KWi. At the end of the device we situate a set of CCD, which measure the intensity of 
light in (y, z) plane. 

Let us describe now how we can resolve, approximatively, any problem 3 without restric- 
tions to the maximum of B + . The key truncation idea can be found in (Ibarra and Kim 
1975), see also (Papadimotriou and Steiglitz 1982). We reduce the case of arbitrary coeffi- 
cients Ci < B + ,Wi to the previously studied case of restricted coefficients. To proceed it, we 
use the binary representations of these numbers. Suppose that each number use < m the 
digits 0, 1, i.e. the size of the binary input is not much than m. Moreover, assume that the 
maximal admissible value of B+ can be written down with only m* digits. Now we trun- 
cate each given number Ci,Wi and B + removing m — m* digits and allowing only first to* 
digits. With these new truncated data we can solve our problem and this solution gives an 
approximative solution. 

The solving procedure to search an approximative solution of the problem 3 can be de- 
scribed as follows. We observe all beams that are measured by CCD and between them we 
choose a beam that have a maximal y-deviation. Such a scheme always gives an approximative 
solution of our problem, with precision of order k/2. 

Let us compare now a deterministic processor and QOD. Let us take the same parameters 
as above. The time processing for our device will be chosen the same as above, i.e., T q = 
W~ 6 sec. For the deterministic computer, according to subsection 2c, one has 

T c ~n 4 Ve-\ 
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where e is the relative precision. The relative error of measurement we can estimate as 
e ~ (k/2)/Rm, which corresponds to the separation of the gaussian beams with amplitudes 
differs by factor like 2. Substituting all the values in the last relation, we have T c 
that much more than T q . So, the quantum optical device gains an acceleration 10 6 times vs 
a typical deterministic computer. 

7 Conclusion 

We have introduced and designed analogous quantum optical devices for simulating dynamic 
programming which provides the following complexity bounds (see notations in subsections 
5.1, 6.1) for different versions of the knapsack problem (section 2). 

Proposition 7.1 i) For the versions 1,2 the implementation cost CI quan t = O(Kn), the 
running time Time quan t = 0{(n + K)M), and the energy cost CE quant = 0{Kn(n + K)M); 

ii) For the version 3 Cl quan t = 0(K 2 n), the running time Time quan t = 0((n + K)M), 
and the energy cost CE quant = 0(K 2 n(n + K)M); 

in) for e- approximative solution of the version 3 CI quan t = 0(n(5 p /e) 2 ), the running time 
Time quant = 0((n + 5 p /e)M), and the energy cost CE quant = 0(n(n + 5/e)(5 p /e) 2 M); (see 
subsections 5.1 and 6.1) 

Also we compare these bounds with ones for deterministic machines (see subsections 5.1, 
6.1). 

Let us discuss briefly the "quantum properties" of QOD. The genuine quantum machine 
has to exploit two quantum properties: i) "exponential resource" connected with exponentially 
large dimension of the state space, and ii ) the "quantum parallelism" which is simultaneous 
evolution in all subspaces of the state space. 

The QOD described above contains only one genuine quantum component, namely, 50% 
mirrors, which split the laser beams. This splitting gives the possibility to realize the "quan- 
tum parallelism". The splitting mirrors operate with the laser beams that are macroscopic 
objects. However, this macroscopity prevents to the realization of "exponential resource". In 
this framework, realization of the exponential resource means using of the mirrors of expo- 
nential size in order to separate the exponentially large number of the laser beams. Moreover, 
in this case we need exponentially large energetic resource in order to keep the intensity of 
the laser beams on a suitable level. So, it is difficult to realize of the "exponential resource 
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with help of macroscopic quantum devices. But the second quantum property - "quantum 
parallelism" - can be realized by described above devices. 
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